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Abstract 

This work reviews foundations and applications of the complex-energy continuum shell model that pro- 
vides a consistent many-body description of bound states, resonances, and scattering states. The model can 
be considered a quasi-stationary open quantum system extension of the standard configuration interaction 
approach for well-bound (closed) systems. 
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1 Introduction 



Small quantum systems, whose properties are profoundly affected by environment, i.e., continuum of scattering 
and decay channels, are intensely studied in various fields of physics (nuclear physics, atomic and molecular 
physics, nanoscience, quantum optics, etc.). These different open quantum systems (OQS), in spite of their 
specific features, have generic properties, which are common to all weakly bound/unbound systems close to the 
threshold. While many of these phenomena have been originally studied in nuclear reactions, it is not possible 
to experimentally control the behavior of the nucleus by varying external parameters as in, e.g., atoms and 
molecules, quantum dots, or microwave resonators. 

Nuclear physics contains the core of sub-atomic science with the main focus on self-organization and stability 
of nucleonic matter. Nuclei themselves are prototypical mesoscopic OQSs and splendid laboratories of many- 
body physics. While the number of degrees of freedom in heavy nuclei is large, it is still very small compared 
to the number of electrons in a solid or atoms in a mole of gas. Nevertheless, nuclei exhibit behaviors that 
are emergent in nature and present in other complex systems. Since nuclear properties are profoundly affected 
by environment, i.e., the many-body continuum representing scattering and decay channels, a simultaneous 
understanding of the structural and reaction aspects is at the very heart of understanding short-lived nucleonic 
matter. An essential part of the motion of those exotic systems is in classically forbidden regions, and their 
properties are profoundly impacted by both the continuum and many-body correlations (see Fig.[T]and Ref. [T|). 
By studying the limits of nuclear existence, we also improve our understanding of the ordinary nuclei around 
us, extending the nuclear paradigm. 
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Figure 1: Schematic diagram illustrating various aspects of physics important in neutron-rich nuclei. One- 
neutron separation energies Sn, relative to the one-neutron drip-line limit (5„=0), are shown for some isotopic 
chain as a function of N. Weakly bound nuclei, such as halos, inhabit a drip-line region characterized by very 
small values of Sn- The unbound nuclei beyond the one-neutron drip line are resonances. Their widths (repre- 
sented by a dark area) vary depending on their excitation energy and angular momentum. At low excitation 
energies, well-bound nuclei can be considered as closed quantum systems. Weakly bound and unbound nuclei 
are open quantum systems that are strongly coupled to the scattering environment. The regime of particularly 
strong many-body correlations involving weakly bound and unbound neutrons with |5„|<2MeV is indicated 
as "correlation dominated" [T]. In this region, the neutron chemical potential has the same magnitude as the 
neutron pairing gap. The astrophysical r-process is expected to proceed in the region of low separation ener- 
gies, around 2-4 MeV. In this region, effective interactions are strongly affected by isospin, and the many-body 
correlations and continuum effects are essential. 
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Resonances are commonly found in various quantum systems, independently of their building blocks and 
the kinematic regime of their appearance. Resonances are genuine intrinsic properties of quantum systems, 
associated with their natural frequencies, and describing preferential decays of unbound states. The effect of 
resonances and the non-resonant scattering states can be considered in the OQS extension of the shell model 
(SM), the so-called continuum shell model (CSM) [2]. A particular realization of the CSM is the complex-energy 
CSM based on the Berggren ensemble, the Gamow Shell Model (GSM). The standard quantum mechanics 
(the Hilbert space formulation) does not allow the description of state vectors with exponential growth and 
exponential decay, such as resonance states. Since the spectrum of an observable is real in the Hilbert space, 
the usual procedure for treating unbound resonance states is either to extract the trace of resonances from the 
real-energy continuum level density or to describe the resonances by joining the bound state solution in the 
interior region with an asymptotic solution, e.g., within the R-matrix approach [3111]. 

These difficulties of the Hilbert space formulation have been resolved in the Rigged Hilbert Space (RHS) 
[U [6] formulation. Thus, in a broader sense, the mathematical setting of GSM follows directly from the 
formulation of quantum mechanics in the RHS (Gel'fand triple) [5l [6] rather than the usual Hilbert space (see, 
e.g., Refs. [7, 8j). The Gel'fand triple framework allows not only to formulate rigorously Dirac's formalism of 
'bras' and 'kets' [9l |10] but also encompasses new concepts, like Gamow states, and is suitable for extending 
the domain of quantum mechanics into the time-asymmetric processes like decays, offering a unified treatment 
of bound, resonance, and scattering states. 

In the past, Gamow states have been discussed in various contexts [TTl[T2l[T3l[M[T5l[T6l[T7l[T8l[T9l[20l[2Tl 
[22] . including those related to the GSM works discussed in this paper. The subject of this review is to present 
the foundations and selected applications of the nuclear SM in the Berggren ensemble. This model provides a 
natural generalization of the standard nuclear SM for the description of configuration mixing in weakly bound 
states and resonances. 

Alternative Hilbert-space-based description of the interplay between many-body scattering states, reso- 
nances, and bound states can be formulated using the concept of projected subspaces [231 [23]. Realistic studies 
within the real-energy CSM have been presented recently in the framework of the Shell Model Embedded in the 
Continuum (SMEC) |25[I26]. A phenomenological CSM with an approximate treatment of continuum couplings 
has been proposed in Ref. [2^. These models have the advantage of being able to provide reaction observables 
such as radiative capture and elastic/inelastic cross sections [25]. However, they rely on a somehow artificial 
separation of the Hilbert space into bound/resonant and non-resonant scattering parts. This is contrary to 
GSM, in which all states are treated on the same footing. The use of projected subspaces in SMEC (CSM) 
complicates matters significantly when several particles in the non-resonant scattering continuum are consid- 
ered. It is only recently that SMEC has been extended to treat the two-particle continuum [26j in the context 
of two-proton radioactivity |28j. While numerically demanding, these first applications were still subject to 
several approximations, such as sequential or two-body cluster decays. 

The paper is organized as follows. Section [2] discusses the concept of Gamow states and the Berggren basis. 
The complex scaling method, employed in GSM to find Gamow states and regularize integrals, is outlined 
in Sec. [3l Section [H describes properties of the one-body Fock space of the GSM. The many-body GSM is 
described in Sec. [H and examples of applications are given in Sec. [H Finally, Sec. [7| outlines perspectives. 



2 Gamow states and the Berggren ensemble 

The Gamow states |29l [30] (sometimes called Siegert [H] or resonant states) were introduced for the first time 
in 1928 by George Gamow to describe a decay. Gamow introduced complex-energy eigenstates 

En = En- i Y (1) 



in order to explain the phenomenon of particle emission in a quasi-stationary formalism [31]. Indeed, if one 
looks at the temporal part of a decaying state, which is e^i?oi/?ie-rt/(2n)^ g^e notices that the squared modulus 
of the wave function has the time-dependence oc e~'"* , and one can identify F with the decay width that defines 
the half-life of the state: 

^ ?iln2 

Ti/2 = -f^. (2) 
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The resonant states can be identified with the poles of the scattering matrix S{kn), 

in the complex-momentum plane. The bound states are thus resonant states lying on the imaginary-A; axis. 

In 1968, Berggren proposed [32] a completeness relation for single-particle (s.p.) resonant states involving 
a complex-energy scattering continuum: 

^ Un{En,r)un{En,r) + [ dEu{E , r)u{E, r ) = 6{r - /), (4) 

where 

Un{En,r) OiiKr) e"'-' (5) 

and kn = inn {i^n > 0) for bound states and kn = Jn — ii^n {i^n-,ln > 0) for decaying resonances in the fourth 
quadrant of the complex-/? plane. 

In practical applications, it is more convenient to write ^ in momentum space: 

V \Un){Un\+ [ \uik)){uik)\dk = l. (6) 

As seen in Fig. [21 bound and quasi-bound states (e.g, resonances) enter the completeness relation ([6]) on the 
same footing. The completeness relation introduced by Berggren reduces to the traditional one involving bound 
and scattering states if the contour lies on the real /c-axis. 

The radial wave functions corresponding to the 5-matrix poles with Im{k) < diverge as r — > oo; hence, 
the scalar product among basis states has to be generalized. Firstly, a biorthogonal basis is used for the radial 
wave function, i.e. a different basis set for the states in bra and ket positions. Second, radial integrals have to 
be regularized. These extensions stem from the fact that continuum states belong to RHS. The RHS metric 
becomes equivalent to the ordinary metric only for bound states. 

While many important developments took place prior to Berggren's work, e.g., analytic properties of the 
Green's function and the scattering matrix were recognized [13| , and Zel'dovich [33j and Hokkyo [M] proposed 
regularization methods to normalize Gamow vectors; however, the bound and resonant states were thought to be 
non-orthogonal and this caused conceptual difficulties. In Ref. [32], Berggren applied the regularization method 
of Zel'dovich (with a Gaussian convergence factor) to normalize Gamow states and proved the completeness 
relation Another type of regularization was introduced by Romo [35] who used the analytic continuation of 
the norm integral from the upper k halfplane to the resonant state in the lower k halfplane. As an application, 
he solved a two-channel problem involving a Gamow state. Zimanyi introduced yet another way of normalizing 
Gamow states by changing the strength of the generating potential [36^ [37] and applying this technique to the 
coupled Lane equations for complex-energy isobaric analog resonances. Bang and Zimanyi used a Gamow form 
factor in the description of stripping reactions leading to a final resonant state [38]. In 1972, Gyarmati and 
Vertse |39] proved the equivalence of normalization procedures of Zel'dovich and Romo for neutron resonant 
states. They demonstrated the existence of the norm for a proton resonant state and introduced a new 
regularization procedure. This external complex scaling method, applied in GSM, is described in some detail 
in Sec. K% 



2.1 Antibound states 

Antibound (or virtual) states lie on the negative semi-axis of the imaginary k'. — ^/^n {Kn > 0) (see Fig.[2|). 
They have real and negative energies that are located in the second Riemann sheet of the complex energy plane 
[401 l4H |42| |43] . Asymptotically, the radial wave function of a virtual state grows exponentially: 

Wn{En,r) ~ Oliknr) ~ e''-"- = e^"^ (7) 

As often discussed in the literature, it is difficult to provide a physical interpretation to virtual states. In the 
Hilbert space formulation of quantum mechanics, a virtual state is not considered as a state but as a feature of 
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Figure 2: Location of one-body states in the complex momentum plane. The Berggren completeness relation 
([6]) involves the bound states (b) lying on the imaginary /c-axis, scattering states on the L"*" contour (solid thick 
line), and resonant decaying states (d) in the fourth quarter of the complex-A; plane lying between the real axis 
and L"*". In the general expansion of the resolvent, the scattering states along the U contour and all resonant 
states lying above U are included. The antibound states (a) can be included in the generalized completeness 
relation; in this case the contour L"*" has to be slightly deformed (dashed thick line). 



the system (as the second energy sheet is considered unphysical and inaccessible through direct experiments). 
In the RHS formulation, the virtual state can be interpreted both as a vector in the RHS and as a pole of 
the S-matrix. The latter implies that virtual states can be "seen" only close to the threshold; they manifest 
themselves in an increased localization of low-energy scattering states. Consequently, the presence of a virtual 
state at a sufficiently small energy has an appreciable influence on the scattering length and the low energy 
scattering cross section. Classic examples include the low-energy £=0 nucleon-nucleon scattering characterized 
by a large and negative scattering length [42 i [H] . Related to this is an increased localization of real-energy 
scattering states just above threshold [i5] . 

It is instructive to study the energies of antibound states as a function of potential parameters. A detailed 
study of this problem was done by Nussenzveig |46j (see also Ref. [47] for a discussion of the one-dimensional 
case). To this end, it is convenient to fix the geometry (radius and diffuseness) of the potential and follow the 
trajectories of the bound and antibound poles as a function of the well depth Vq > 0. If there is no potential 
barrier, e.g., for s-wave neutrons, no narrow resonances appear. Here, by increasing Vq from the minimal value 
at which the antibound state appears, the pole goes across /c = 0, forming a bound state with a radial wave 
function having one node, n = 1. For protons or for I > neutrons, a potential barrier is present which can 
support resonant poles. In this case, decaying and capturing poles appear in pairs which move towards the 
imaginary /c-axis as Vq increases. At certain values of Vq the poles meet at k=0, forming a double singularity. 
For still larger values of Vq, one of the poles becomes a bound state, while the other one moves down as an 
antibound state. 

Neutron /=0 pole trajectories are shown in Fig. [3] as a function of Vq for a square well potential. The 
behavior of the antibound poles conforms to the general pattern discussed above. Namely, below Vq = 1.05 
MeV, there appears only one antibound state. This pole moves to the upper half fc-plane and forms a bound 
state with one node. By increasing the depth further from Vq = 9.05 MeV, there appear two antibound poles. 
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Figure 3: Positions of the bound and antibound neutron /=0 poles with n=l,2,3 radial nodes on the imaginary 
/c-axis as a function of the depth Vq of the square well potential with radius R = 7 fm. 



One of them moves upwards and at about Vq = 9.5 MeV, it emerges as a 2s bound state, while the other 
one moves down and remains antibound. At about Vq = 26 MeV three antibound and two bound states are 
present. By increasing Vq to 30 MeV, the antibound pole lying closest to the origin becomes bound, resulting 
in three bound and two antibound states. At Vo=30 MeV, the 3s bound state is closer to the origin than the 
corresponding antibound pole. Therefore, there are no antibound poles between the loosest bound state and 
the threshold [IQ]. For the cases with a non-zero barrier, the resonant poles meet in the origin and one obtains 
a similar pattern as in Fig. ([3]) but with the Im(/c)-axis shifted down by about A; = 0.15 fm~^. This shift would 
change the relative distance of the antibound and bound states with respect to the origin, with one antibound 
pole to appear between the loosest bound state and the threshold [30] • 

2.2 Generalized Berggren representation 

In the original formulation of the completeness relation ^ , only decaying poles with 

arg{kn) > Re{En) > (8) 

were considered. Consequently, virtual resonances (resonant poles with Re{En) < 0) and antibound states 
were excluded from the Berggren ensemble. The Berggren representation has been generalized in Ref. [IS] by 
including a virtual state into the discrete part of the basis. The generalization requires a slight deformation 
of as indicated in Fig. [2j More recently, this generalized Berggren representation was used in a full form 
(i.e., including the scattering component) for the description of drip-line nuclei ^^Li and ^^Ca [491 I50| . A 
large destructive interference between the virtual state and the scattering contour states has been observed. 
In this picture, the bound ground state of ^^Li could be explained in terms of unbound states of ^''Li. Soon 
afterwards it was noted ^51j that a generalized Berggren basis is not optimal as far as the number of basis 
states is concerned. Indeed, as the halo wave function has a decaying character at large distances due to the 
exponentially increasing asymptotics of the virtual state ([7]), its inclusion in the basis always induces strong 
negative interference with the scattering states. Therefore, adding antibound states to the basis is not beneficial 
in GSM applications, as more discretized scattering states are necessary to reach a required precision. 
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2.3 Berggren completeness relation for protons 

The completeness relation by Berggren ([6|) can be derived from the Newton completeness relation for the set 
of real-energy eigenstates of a short-range potential [40]: 

r+oo 

V|u„)(u„|+/ \uik)){uik)\dk = l. (9) 

The proof can be carried out by deforming the real momentum axis associated with positive-energy scattering 
states into the L"*" contour of Fig. [2] and applying the Cauchy integral theorem [11]. This is possible because 
resonant states appear as residues, due to their S'-matrix pole character. 

Bound and real-energy scattering states generated by a potential bearing a Coulomb tail naturally form 
a complete set as well, as all self-adjoint operators possess a spectral decomposition [52]. However, methods 
used in Ref. [52], relying on the Lebesgue measure theory, are very general and fairly abstract. Thus, a simple 
demonstration of the completeness relation for potentials with a pure Coulomb asymptotic behavior is called 
for. The proof of Q relies on the analyticity of the Green's function of a local potential that is integrable for 
r — > +00 i40|. Potentials having a Coulomb tail, proportional to r~^, have not been considered in Ref. |40j . 
Moreover, non-locality often appears in practical applications [53], so it is important to demonstrate the 
completeness relation for non-local potentials having a Coulomb asymptotic behavior. This has been done 
recently in Refs. [531 [S^, where detailed derivations can be found. 

In short, the main assumption behind the proof presented in Refs. |531 154j is that beyond a finite radius i?o 
the potential behaves as r~^. The corresponding radial Schrodinger equation reads: 

u"{k, r) = (^-^^^^ + v{r) - k"^^ u{k, r) + w{r, r') u{k, r') dr\ (10) 

where v{r) and w{r,r') are respectively the local and non-local potentials. For r > Rq, the wave function 
is a linear combination of regular and irregular Coulomb functions. Their analytical character and the fact 
that regular Coulomb wave functions form a complete set [Ml [S51 ES] are important to the demonstration of 
the completeness. Using analytic continuation arguments, one can also demonstrate that the set of bound 
and real-energy scattering states generated by complex potentials v{r) and w{r,r') is complete, provided no 
exceptional points, i.e., bound states of norm zero [57j, are present. 

2.4 Foundation of the GSM and interpretation of resonant states 

Heuristic methods often precede a complete formulation of the physical theory in terms of an adequate math- 
ematical apparatus. This was the case for Dirac's formulation of quantum mechanics [58] which a posteriori 
found a satisfactory setting in the RHS [9l [1^. This is also the case for theoretical developments utilizing 
Gamow states to describe weakly bound and/or unbound states of quantum systems. The Berggren com- 
pleteness relation [32], which replaces the real-energy scattering states by the resonance contribution and a 
background of complex-energy continuum states, puts the resonance part of the spectrum on the same footing 
as the bound and scattering spectrum. However, resonances do not belong to the usual Hilbert space, so the 
mathematical apparatus of quantum mechanics in Hilbert space is inappropriate and cannot encompass con- 
cepts such as Gamow states. Quite unexpectedly, it turned out that the mathematical structure of the RHS is 
ready to extend the domain of quantum mechanics into the time-asymmetric processes such as decays. 
The RHS, or Gel'fand triple, is a triad of spaces [5l iGj [59]: 

^cnc^'', (11) 

which represent different completions of the same infinitely dimensional linear space ^. <I> (the subspace of 
test functions) is a dense subspace of Hilbert space Ti., and (the space of distributions) is the space of 
antilinear functionals over $ is also the largest subspace of TC on which expectation values, uncertainties 
and commutation relations can be correctly defined for unbounded operators. 

The linear functionals over $ are contained in the space <I> which is related to another RHS: 

$ C W C . (12) 
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One should stress that 7i in (jlip and (jl2p does not have any particular significance in the RHS formalism. 
For all physics problems, one needs only the dual pair of spaces $ C which characterize the considered 
quantum system. For example, bras and kets, which are related to the continuous part of the spectrum 
of an observable, belong to <1> and <^^, respectively, and not to TL. Apparently, the RHS provides also a 
more convenient framework than 7i to capture physical principles of quantum mechanics. The Heisenberg's 
uncertainty relations are properly defined on and not on Ti. 

Besides implementation of Dirac's formalism of bras and kets, the RHS provides a framework for a quantum 
mechanical description of common irreversible processes, like the formation or decay of quantum states. Indeed, 
$^ may contain generalized eigenvectors of the observable(s) with complex eigenvalues. The generalized eigen- 
vectors of the Hermitian Hamiltonian are Gamow vectors. Examples of the RHS for the observables having a 
continuous part in the spectrum can be found, for example, in Refs. I22j. 

Gamow vectors are state vectors of resonances. They belong to RHS and not to Ti since the self-adjoint 
operators in 7i can only have real eigenvalues. Like the plane waves (the Dirac kets), the Gamow functions are 
not square integrable and must be treated as distributions. In this way, one succeeds to generate Gamow bras 
and kets which fit naturally in the RHS. However, unlike the Dirac kets, the probability density for Gamow 
functions is not constant but increases exponentially in space. 

Rules for the normalization of Gamow functions and corresponding completeness relations have been pro- 
posed [Ml EH EHl EH [62l [3^ independently of the RHS formulation. Similarly, the momentum space represen- 
tation for Gamow functions has been introduced [631 [M] . RHS offers a mathematical framework to formulate 
those results rigorously. In particular, the RHS provides a unifying dual description of bound, resonant and 
scattering states both in terms of the S-matrix and in terms of the vector in the RHS. The transition amplitude 
A{En E) from a resonance of energy En to a scattering state at an energy E {E > 0)[22j: 

A{En -^E)= iV2^Nn5{E - En), (13) 

where Nn denotes the normalization factor, is proportional to the complex (5-function [65] which, far off the 
threshold energy in the region of the resonance, can be approximated by the Breit-Wigner amplitude. This 
physically expected result provides a formal link between Gamow states and nearly-Lorentzian peaks in cross- 
sections for narrow, isolated resonances. 

Interpretation of certain aspects of the RHS formulation of quantum mechanics is still debated, and con- 
nections between Berggren and RHS formulations continue to be studied (for a recent review, see [60]). One 
aspect of this discussion concerns complex matrix elements of the operators and probabilistic interpretation of 
the resonant wave function. As an illustration of this problem for resonances, let us consider the solution of 
the Schrodinger equation x(r, t) for a s.p. Hamiltonian h: 

in^^x{r,t) =hx{r,t). (14) 

In the stationary picture, the resonant wave function is the product of time-dependent r(t) and space-dependent 
tp{r) factors: 

X(r,i) =r(t)V'(r), (15) 

where the time-dependent factor is: 

r(t) = e-'^K (16) 

The complex energy En ([I]) is the generalized eigenvalue of h with the complex wave number kn ([3]). The 
corresponding space-dependent eigenfunction ip{r,kn), 

= ^niM-^ kn) = ^^^^^^^^m)XsU, (17) 

satisfies an outgoing wave boundary condition. At large distances, the radial part of the i=0 wave function is 
proportional to exp{iknr). For a decaying resonance, kn = Jn — it^n {in, > 0), the time-dependent factor 

T(t) = e-'—* = e"'— *e- , (18) 
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describes an exponential decay with a decay width r„,. As r increases, the radial part of ip{r, kn) oscillates with 
an exponentially growing amplitude: 



The decay process in this picture is unlimited in both time {t £ [0, oo]) and space (r G [0, oo]) and, henceforth, 
continues forever. This feature guarantees, in fact, the particle-number conservation [66]: the exponential 
temporal decrease of the wave function amplitude is complemented by its exponential spatial increase, and the 
divergence of the resonance wave function assures that the particle number is conserved. 

In general, the quasi-stationary description of a time-dependent physical process works well when the 
decaying state is narrow, as in this case the formation of a state and its subsequent decay can be well separated 
in time. A narrow Gamow resonance has a large overlap with a properly shaped wave packet [67] and its radial 
wave function is similar to a bound state wave function while the imaginary part is small inside the nucleus. 
In the reaction cross section, such a narrow resonance shows up as a sharp peak which can be well separated 
from the non-resonant background. 

Another aspect of the debate on the physical interpretation of Gamow vectors concerns the expectation 
value of an operator in a resonant state. Berggren's interpretation of real and imaginary parts of the expectation 
value of an operator is based on the discussion of interference effects in the reaction cross-sections |68| . Berggren 
showed that the imaginary part of the complex cross section calculated with the Gamow resonance as a final 
state describes the interference of that resonance with the non-resonant background. In a later paper [69], 
Berggren demonstrated that for any operator which commutes with the Hamiltonian, the expectation value of 
that operator is the real part of its matrix element, and the negative of the square of the imaginary part of the 
matrix element can be associated with the square of an uncertainty of the expectation value. This proof holds 
only for operators which commute with the Hamiltonian. It turns out that the Berggren interpretation of the 
expectation value [69j is largely equivalent to the RHS formulation by Bohm and Gadella |70| and coincides 
with it in the leading order [151 E] • 

The meaning of a complex radial expectation value was studied in Ref. [72] using a model of two coupled 
channels [73]. It was found that the radial expectation value in this model can be interpreted according to 
Berggren's prescription and that it can be used to characterize the resonances as long as the wavelength of the 
decaying wave is shorter than the radial extension of the resonant state. 

2.5 Applications of single-particle resonant states to various nuclear structure problems 

Already in 1970, in the pioneering work employing resonant states, a truncated resonant basis composed of 
bound neutron and resonant proton states was used to describe isobaric analog resonances [37]. The pole 
expansion of the radial wave function, the Green's function, and the scattering matrix were subsequently 
discussed in a rigorous way in Ref. [74j by using the Mittag-Leffler (ML) expansion. They used the normalization 
of the pole solutions introduced in [3l] and later employed in, e.g., Refs. [611 [75]. The authors of Ref. [71] 
suggested that the pole expansion could be used in (shell-model-type) nuclear structure calculations; however, 
they worried about the non-orthogonality of the basis. In a numerical study [76], the usefulness of the ML-type 
pole expansion of the Green's function was examined. The question of overcompleteness was finally put to rest 
by Berggren and Lind [77] . 

Resonant state expansions were also used by atomic physicists for calculating the S'-matrix, expressed as a 
sum of pole terms and a contour integral [78]. The positions of the poles were calculated by using the uniform 
complex scaling (see Sec. [3|), while the residues were obtained by using numerical methods described in Refs. 
[791 [80] . Neutron emission from a heavy-ion reaction was modeled in a time-dependent two-center shell model 
in Ref. [81]. Therein, resonant states were calculated in the momentum space using the potential separable 
expansion (PSE) method [821 183]. in which the potential is expanded in a square integrable (harmonic oscillator; 
HO) basis, and the proper asymptotic of the wave function is guaranteed by the presence of the free Green's 
function. 

A less severe truncation of the Berggren basis was employed in the late 1980s in the resonant approach to 
the random phase approximation (RRPA), in which the contribution from the non-resonant continuum was 
neglected [HH |85]. Such an expansion is usually referred to as a pole approximation. (A detailed summary 




ik, 
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of the early theoretical approaches employing resonant states can be found in Ref. [86j.) A similar pole 
approximation was used in the RRPA description of giant resonance escape widths ^87]. A careful test of 
different pole expansions was carried out in Ref. [76J by using a benchmark continuum RPA result for the 
particle-hole response function (see Sec. 12.61 for details). It was found that the ML-type pole expansion was 
closer to the continuum RPA, based on the Green's function, than the pole approximation. This result was 
explained by Berggren and Lind [77j who identified the ML-type pole expansion of the resolvent by using the 
so-called U contour (see Fig. [2]) in the complex /c-plane. Unfortunately, an expansion using the U contour, 



cannot be used in shell-model-like calculations because the factor 1/2 in the front of the sum over the resonant 
(bound, antibound, capturing, and decaying) states destroys the idempotency of the unity operator [88j. In 
addition, pole expansions of type (j20l) are accurate only if very many poles (up to some large energy cut-off) 
are included in the sum; contributions from broad and virtual resonances (having negative real energies) are in 
fact essential [88] . 

In Ref. [89], Gamow states were expanded in a finite r £ [0,Ro] interval by using a finite number of 
square integrable basis states, called Siegert pseudostates (SPS). The expansion coefficients were determined 
by solving a normal eigenvalue problem where the dimension of the system has been doubled. This doubled 
system represents an eigenvalue problem with a weight matrix of the Bloch operator. The weight matrix allows 
the continuation of the SPS into the tail region and defines a special inner product in which a completeness 
relation can be constructed for the SPS. A pole expansion of the Green's function, similar to the ML expansion 
of Ref. [88], was carried out. This expansion was further exploited in Ref. [90]. The usefulness of SPS was 
further demonstrated in studies of time evolution of an infinitely extended system [91] . 

The resonant and continuum components of the strength function were calculated by using the Berggren 
completeness relation in studies of Coulomb breakup of ^^Be |92) . They found the non-resonant contribution 
of the part of the contour which returns to the real fc-axis negligible. A similar result was obtained in Ref. [93] . 

In some of the pole expansions mentioned above, complex eigenvalues and radial wave functions of resonant 
states were calculated using the code GAMOW [9l]. Later, more efficient numerical techniques have been used 
based on the piecewise perturbation method [951 [96] . In other work [97] spherical resonant states were obtained 
by solving a non- linear set of equations. Bound and resonant state energies were calculated in code pseudo 
by using the PSE method in [98] for an axially symmetric potential by finding the zeroes of the Fredholm 
determinant of the Lippmann-Schwinger equation. Recently a method was developed in which the momentum- 
space Schrodinger equation was solved by a contour deformation method [93] . 

More recent applications of pole expansions include a description of giant multipole resonances [99] (see 
Sec. 12. 6p . and partial decay widths corresponding to the proton decay from the Gamow- Teller and isobaric 
analog state resonances [lOOj . The Berggren representation in its full form, including the complex scattering 
continuum, was employed in Ref. [lOlj to describe the eigenstates of a realistic one-body nuclear Hamiltonian. 
In this work, eigenstates of a spherical symmetric potential were expanded in a full Berggren basis of a different 
potential. (Examples of such expansions are presented in Sec. 14.51 ) The full Berggren ensemble was employed 
in studies of the average s.p. level density [102] and shell corrections [1031 110 j] . Quasiparticle resonances can 
be constructed from the Gamow resonances in the BCS approximation by generalizing the standard BCS to 
the complex energy plane if the contour integral is neglected. If the non-resonant contribution is included, the 
contribution from the continuum level density has to be added [105] when calculating the particle number. 
However, a proper description of pairing correlations in weakly bound systems must go beyond the BCS 
approximation [106| . Here, the tool of choice is the Gamow-HFB method to be mentioned in Sec 14.41 

Gamow states and resonant expansions have been extensively used in the description of alpha, proton, 
and neutron emission from nuclei. For very narrow alpha and proton resonances, a single-channel Gamow 
approximation that ignores the non-resonant background is fully adequate [107[ I108j . An attempt was made 
to calculate resonant states of a deformed potential by means of an expansion in a spherical Berggren basis 
[10911110] . However, for a description of narrow proton resonances in deformed nuclei, the accuracy of spherical 
expansion turned out to be not satisfactory and alternative approaches were used. One of them was the direct 
numerical integration of the set of coupled differential equations using the piecewise perturbation method. 
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The resulting code CCGAMOW was employed in Ref. [Ill] to explain full and partial decay widths of protons 
from axially symmetric nuclei within an adiabatic approach that ignores the rotational excitations of the 
daughter system. A full non-adiabatic description of proton emitters was carried out in Refs. [112\ 1113] using a 
modified version of the coupled channels program nonadi, and the extension to the triaxial case can be found in 
Ref. |114j . A shell-model analysis of the proton emission from ^^Cl using Gamow wave functions was performed 
in Ref. [115j . In a study of Ref. |116j . Gamow states were used to benchmark a modified two-potential approach 
to tunneling problems. 

The applications of resonant states mentioned so far were concerned with the s.p. space. By 2003, however, it 
became possible to apply the Berggren ensemble to a general many-body case. In two parallel studies [11711118] . 
based on a configuration mixing approach employing the one-body Berggren ensemble, two-particle resonances 
in ^°Ni and ^'^^Te were calculated using separable multipole interactions [117] . and systematic calculations for 
the neutron-rich nuclei ^"^'^He and 18-22q -^gj-g carried out using the surface-delta interaction [119] . More 
applications soon followed (see Sec. [6]). In the following, we shall refer to these two realizations of the SM in 
the complex energy plane as the Gamow Shell Model since the underlying concept is the same. 



2.6 Example 1: resonant state expansions and continuum RPA 

Continuum RPA (CRPA) is a useful tool for testing the usefulness and accuracy of resonant state expansions. 
For separable multipole-multipole interactions, the solution of the CRPA equations is obtained by finding 
complex roots of the corresponding dispersion relation \76\ 199] . Complex energy poles of the particle-hole 
response R{E) yield positions and full widths of the correlated resonances while the partial decay widths of 
these states can be obtained from the associated residua. 

Within CRPA, the particle Green's function in a selected partial wave is 

/ / ,x -u(fc,r<) v{k,r^) 

^(^'^'^) = w{^) — ' 

where W denotes the Wronskian of the regular u{k,r) and irregular v{k,r) solutions at r<=min(r, r') and 
r>=max(r, r'), respectively, at real values of k. In the standard Berggren expansion using the L contour of 
Fig. [21 the corresponding Green's function can be written as: 

/ / , N Uiiki,r) Ui{ki,r') 1 f u{q,r)uiq,r') 

gL{r,r ,k) = "2 72 h - / T- r — dq. (22) 

if^d k^-kf t^Jl q{k-q) 

On the other hand, if one uses the U contour of Fig. [2l the Green's function takes the form: 

. / ,^ Ui{ki,r) Ui{ki,r') 1 f u{q,r)u{q,r') 



The two complex-energy representations of Green's function given by Eqs. (j22|23p and the real-energy expression 
(|2ip are mathematically equivalent; hence, they should give the same results in practical applications. 

Prom a practical point of view, it is interesting to know if one can choose the complex contours in such a 
way that the contributions from the integrals in Eqs. (j22p and (|23p become negligible. This has been checked in 
Ref. ^76j , in which the particle-hole response functions were calculated for a square- well potential supplemented 
by the Coulomb field. For narrow and isolated resonances, the partial widths obtained by using the leading pole 
term in Eqs. (|22]l and (f23]l agree with the exact result within 10 %. The subsequent calculations of Refs. [881199] 
using the Woods-Saxon potential demonstrated that the contribution from the U contour to the p-h response 
function could be further reduced but this does not hold for the L contour. Furthermore, it was shown [771 [H] 
that neglecting the integral in Eq. (|22p destroys the symmetry of the resolvent and causes an artificial threshold 
behavior of the Green's function. 



2.7 Example 2: GSM for isobaric analogue resonances in Lane model 

The isobaric analogue resonance (lAR) can be described phenomenologically by the coupled-channel Lane 
equations (CCLE) The lAR-CCLE problem offers an excellent opportunity to compare the complex-energy 
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GSM approach with the standard CCLE sokition along the real energy axis. Such a test has been carried out 
in Ref. [120J for the Lane Hamiltonian. Within CCLE, the multichannel Schrodinger equation has been solved 
by a direct numerical integration. The resulting proton scattering S-matrix has been calculated along the real 
energy axis in the region of the lAR and then fitted by using a single pole approximation, 



Ep - SxATl , 



(24) 



to determine the position and width F of the lAR. In Eq. (j24|) . the background phase shift 6p{Ep) was 
assumed to have a linear energy dependence in order to better reproduce the non-resonant background. The 
best fit parameter values, denoted by E^{CCLE) and T{CCLE), are listed in Table[TJ They coincide with the 



GSM eigenvalue Sjar = Er{GSM) 



■ T(GSM) 



within 1 keV even for the broad states. The integrated effect of 



Ij 


Er{GSM) Er{CCLE) 


T{GSM) T{CCLE) 


Sl/2 

d^/2 


16.445 16.444 
16.918 16.917 
17.441 17.440 


0.141 0.140 
0.156 0.156 
0.144 0.145 



Table 1: Comparison of the lAR parameters calculated by using CCLE and GSM for three partial waves. All 
energies are in MeV. 

the proton continuum along the complex path is small but essential to yield the correct value of T{GSM). 

3 Complex scaling method 

Another approach that employs the concept of complex energies and distorted contour is the complex scaling 
method (CS). Some aspects of the CS formalism, e.g., the use of complex coordinates, are very relevant to the 
GSM. For a comprehensive review of the CS method, we refer the reader to Reinhardt |121j . Below only the 
most essential facts are summarized. 

Within CS, like in GSM, the resonances are represented by poles of the scattering matrix S{k). By making 
use of the CS transformation on particle coordinates, 

r ^ e^V, (25) 

one guarantees that wave functions of the selected resonances become square integrable |122l 11231 I124j . A 
unitary CS operator U (9) acting on a s.p. wave function gives: 

f7(0)^(r) =e^iV(re^^), (26) 

where the factor e*2^ comes from the three-dimensional volume element [121] . Under U{0), the Hamiltonian 
transforms as 

heir) = U{e)h{r)lJ{6)-\ (27) 

(This variant of CS is referred to as uniform complex scaling] it is to be distinguished from the exterior complex 
scaling discussed below.) The transformed Hamiltonian hg is no longer hermitian as it acquires a complex 
potential. However, for a wide class of local and nonlocal potentials, called dilation- analytic potentials, the 
so-called ABC theorem [T22 | [T23 | [T23]) is vahd. This theorem states that: 

i. The bound states of h and Kq are the same; 

ii. The positive-energy spectrum of the original Hamiltonian h is rotated down by an angle of 29 into the 
complex-energy plane, exposing a higher Riemann sheet of the resolvent; 

iii. The resonant states of h with eigenvalues En satisfying the condition \arg{En)\ < 29 are also eigenvalues 
of hg and their wave functions are square integrable (see Fig. U]) . 
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Figure 4: The complex energy plane showing bound poles (dots), resonant poles exposed by the complex 
scaling transformation (asterisks), the rotated continuum of the scattering states (dotted line), and a hidden 
pole (marked by a square; not exposed by the CS transformation). According to the ABC theorem, the 
eigenfunction of the exposed poles are square integrable. 

Unfortunately, some of the phenomenological potentials widely used in nuclear physics are not dilatation- 
analytic (e.g., the Coulomb potential of a uniformly charged sphere) or they are dilatation-analytic only in a 
limited range of 9. For instance, the commonly used Woods-Saxon potential and its derivative are dilatation- 
analytic only for 9 < 6c = arctg{^). Some of these difficulties, however, can be cured by the use of the 
exterior complex scaling (ECS), originally introduced in Ref. [39j to prove the existence of the norm of a 
Gamow resonance for charged particles (see also Ref. |125] ). Within ECS, the rotation of r to the complex 
starts at a finite distance at or beyond the range of the nuclear potentials > R. The operator Ua{9) of 
the exterior complex scaling is defined for a wave function ^(r) by 



The ECS influences only the asymptotic behavior of functions and can be applied to any finite range potential, 
including the Coulomb potential in the region where it behaves as r~^. 

The Berggren ensemble used in GSM contains non square-integrable functions; hence, a regularization 
procedure is needed to calculate norms and matrix elements needed in the GSM. The ECS can be used for this 
purpose. The complex-scaled states have square integrable wave functions and all related integrals are finite. 
Since the resulting matrix elements are independent of the angle, they can be calculated at any convenient value 
of 9. However, if one is interested in the radial dependence of the wave function for real r (e.g., to estimate 
partial widths), the CS transformation needs to be inverted. 



Such a back rotation can introduce large errors |126j if the solution is expanded and approximated by a finite 
number of terms. If Fade approximants are used for performing back rotation |127j . good accuracy can be 
achieved for narrow resonances. A comparison of the diagonalization in a Berggren basis and the CS method 
resulted in excellent agreement for the complex energy of the lAR in the Lane model ^120j . 




(28) 



V; = C/(0)-Ve. 



(29) 
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The CS and ECS methods have been widely used in many-body calculations of unbound states and in the 
description of atomic and nuclear reactions, see, e.g., Refs. [ml Ell \lM IHIl \IM IHSl [iMl \IM [l36]. For 
recent nuclear examples, see also Sec. [7| 



4 One-body space of the GSM 

This Section describes the s.p. space of the GSM, i.e., the resonant and complex-/c scattering states generated by 
an auxiliary s.p. Hamiltonian. We describe the normalization procedure for the continuum states and discuss the 
discretization of the scattering contour. The optimal basis-generating potential is obtained from the Hartree- 
Fock procedure generalized to resonant states. Finally, selected examples of expansions in the Berggren basis 
are given. They nicely demonstrate the completeness of the Berggren ensemble and the accuracy of calculations. 

4.1 Resonant and scattering states of generating potential 

The radial s.p. wave functions u{k, r) that enter the Berggren ensemble dGj) are solutions of the radial Schrodinger 
equation: 

^' ^ ^ +^V{r)-k^]u{k,r), (30) 



where i is the orbital angular momentum, m is the reduced mass of the nucleon, and V{r) is the one-body 
potential that generates the basis. The choice of the potential is in principle arbitrary, but in practice, one is 
trying to optimize it by adjusting its parameters to experimental s.p. states. Another option to calculate it is 
using the Hartree-Fock (HF) method, as discussed in Sec. 14.41 

In many applications, a Woods-Saxon (WS) central field is used, supplemented by the spin-orbit and 
Coulomb potentials: 



V{r) = -Vofir)-4Vsoil-s)- 

r 

where 

-1 



df{r) 



dr 



+ Vcir), (31) 



f^r)={ l1 + «^p("^)]" ,ifr < i? (32) 

, if r > , 

is a WS form-factor, Vq is the WS potential strength, Rq is the radius, d is the WS potential diffuseness, R is 
a cutoff radius, Vso is the spin-orbit coupling strength, and Vc{r) is the Coulomb potential for protons. The 
solution u{k,r) is regular at the origin, i.e., 

u(A;,r) ~ Co/+i , r ^ 0. (33) 

At large distances, where the nuclear part of the potential is zero, u{k, r) satisfies the asymptotic radial equation: 

d'^u{k,r) f£{i+l) 



dr'^ V r"^ 



+ ^-kAu{k,r), (34) 
r ) 



whose solutions are the regular and irregular Coulomb functions Fi{kr,rj) and Gi{kr,rj), respectively, and 

rj = ^^^^ is the Sommerfeld parameter. The outgoing and incoming Coulomb waves are usually denoted as 
kh 

= Gi{kr,rj) it iFi(kr,r]). For neutrons, = and the asymptotic solutions are the Hankel functions. The 
boundary condition aX r > R can be written as 

u{k,r) ~ C+ff^(A;r) , r > i? (resonant states), (35) 
u{k,r) ~ C+i7^(/i;r) + C_ff^^(/cr) , r > i? (scattering states), (36) 

where Co,C+, and C_ are normalization constants. For the resonant states, i.e., the S'-matrix poles, the wave 
functions have a purely outgoing character; they form a discrete component of the Berggren ensemble. 
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In most applications, Gamow states are expressed in the coordinate representation; they are found by direct 
integration of the Schrodinger equation with outgoing boundary conditions. This procedure is straightforward 
for spherical local potentials and can be applied to deformed potentials by using a coupled-channel decomposi- 
tion of the wave function |112j . In some cases, however, it is useful to work with Gamow states in momentum 
representation [63\ I64j. The continuum discretization in /c-space can be carried out using the Fourier-Bessel 
transform. The Schrodinger equation in momentum representation reads: 

u(k)+ V(k,k')u(k') dk' = Eu(k), (37) 

2m Jo 

where V{k, k') is the Fourier-Bessel transform of the generating potential. Contrary to the coordinate repre- 
sentation, where non-local potentials have to be treated with the equivalent potential method |137j . local and 
non-local potentials are treated on the same footing in the momentum representation. Moreover, discretiza- 
tion of the Schrodinger equation can be replaced by matrix diagonalization |138j . Recent applications of the 
momentum representation can be found in Refs. |139i 193]. This way of generating of Gamow states turned out 
to be very effective, especially for well-deformed nuclei. 



4.2 Normalization of single-particle states 

The s.p. states in (pO|) need to be properly normalized, either to unity (resonant states) or to Dirac delta 
(scattering states). 

The normalization to the Dirac delta is straightforward to achieve, as it is equivalent to the condition 
27rC-|_C_ = 1 in Eq. (j36p . which comes from the asymptotic form of Hankel/Coulomb wave functions. The 
values of the constants C_)_ and C_ can be obtained from the matching condition &i r = R: 



d 
dr 



C+H+{kR) + C^H^^ikR)] = -^u{k, R), (38) 



where u{k, R) is determined by the direct integration from r = to r = R. 

The normalization of decaying states is, however, more complicated as, contrary to bound states, they 
diverge in modulus for r +oo. A solution to this problem can be obtained by considering the analytic 
property of the norm of u{k,r), denoted as M{k). If k is positive imaginary, i.e., the state u{k,r) is bound, 
M{k) is given by 



M{k) = J / n(A;,r)2 dr. (39) 





For resonant states with a nonzero real part, the exterior complex scaling described in Sec. [3] can be used: 

rR r+oo 



pn n+oo 

u{k,rf dr = J u{k,rfdr + ClJ H+{krfdr 
rR r+oo 

= J u{k,rfdr + Clj H+^{kR + kxe'^f e'^ dx , (40) 



TT 

where < < — is the angle of complex rotation. Analyticity of u(/c, r) with respect to r implies that the 

integral (j40p is independent of R and 6. Indeed, if u{k, r) is a decaying state, there always exists a minimal 
angle > for which the complex integral of Eq. (j40p converges for 6 > to a value independent of R and 
9. Hence, one can define the analytic continuation of M{k) as the square root of Eq. (|40p. 
The normalized s.p. wave function V'(^) r) becomes: 

^(fc,r) = ^^[y,(r)x4„. (41) 

All the angular momentum/isospin algebra of the GSM is carried out exactly in the same way as in the standard 
SM. The difficulties inherent to the Berggren basis are all contained in the radial wave function u{k,r). 
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4.3 Discretization of the scattering contour 



In practical applications, as the s.p. basis must be finite, the contour integral along L"*" is limited to the finite 
range k < kmax and then discretized according to a preferred quadrature. The momentum cut-off kmax must 
be sufficiently large to provide convergence. The resulting discretized completeness relation reads: 

\un){un\+ \u{k)) {u{k)\dk ^ \Un){Un\ +YuJi\u{ki)){u{ki)\ =Y\ui){ui\, (42) 
ne(fe,d) ne{b,d) i i 

where {ki,uJi) is the set of discretized momenta and associated weights provided by a quadrature (usually, a 
Gauss-Legendre quadrature is chosen). The discretized scattering states 

Ui{r) = ^iu{ki,r) (43) 

are normalized to unity as the discretization implies the replacement of 5{ki — ki') by ooiSai. The relation (j42p 
involves both resonant and discretized scattering states; it is formally identical to the standard completeness 
relation in a discrete basis. 



4.4 Gamow Hartree-Fock potential 

The nature of s.p. resonant states entering the Berggren ensemble changes with the depth of the generating 
potential. Obviously, a broad resonant state becomes progressively narrow and eventually bound as the po- 
tential's depth increases. This implies that, in order to optimize the s.p. basis, the generating potential should 
be chosen according to the nucleus studied. Let us consider an example of the He chain, treated in a ^He 
core-l-valence particle framework, discussed later in this paper. The potential simulating the '^He core is fitted 
to reproduce the s.p. resonances in ^He. While the resulting basis will work well for a two-neutron halo nucleus 
^He, it is not well suited for a four-neutron halo system ^He. 

An optimal generating potential can be obtained by means of the HF method from the underlying many- 
body Hamiltonian. However, in order to be used in the completeness relation Q, the HF potential has to be 
spherical and the HF approach has to be extended to unbound systems. This approximation is referred to as 
the Gamow Hartree-Fock (GHF) [53]. The starting point of GHF is a GSM Hamiltonian 

HGSM = Y.^hi + Y^Vij, (44) 

i i<j 

consisting of a one-body term 

^ = £^ + (45) 

and a residual two-body interaction Vij acting among valence particles. The one-body field U{r) in Eq. (j45p 
represents the core- valence interaction. It can be taken in a WS form (j3ip . In order to impose spherical 
symmetry, a uniform-filling approximation can be employed, in which nucleons occupy uniformly the valence 
shell (i.e., no individual HF orbitals are blocked). The matrix elements of the HF potential Uuj between two 
spherical states a and f3 carrying quantum numbers {£, j) are: 

{a\U^f\P) = {a\h\P) + ^ J^i^{am Xmx\V\Pm Xmx), (46) 

where A is an occupied shell having quantum numbers N(\) is the number of particles occupying the 

A-shell, and m, mx run over all possible values. One recalls that Uuf provides the exact HF potential only for 
closed-shell nuclei, with which all Ip-lh excitations from the GHF ground state vanish. 

If the nucleus has only one particle or hole outside closed shells, it is possible to define a better approximation 
to Uuf, which has been called M-potential in Ref. [53] and is denoted as Um- In this case, one blocks the last 
particle or hole in the one-body state of largest angular projection m, so that the HF ground state is coupled to 
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J=m. The spherical M-potential is defined by averaging the HF potential over the magnetic quantum number 
m: 

1 ^ 

{apMW) = {a\h\P) + — E AmA), (47) 

m=j+l-Nij X,mx 

where Nij is the number of nucleons occupying the valence shell It is straightforward to show that all 

Ip-lh excitations from the GHF ground state involving the particle in the unfilled shell vanish in this case. 
While this does not hold for the hole states, one can still expect the M-potential to perform better than Uuf 
in this case. 

If the HF ground state is unbound, one is forced to disregard the imaginary part of the GHF potential. For 
the case of one valence particle, the imaginary component of the GHF potential disappears due to antisym- 
metrization, and the imaginary part generated by the core shells must be spurious, as the core is assumed to be 
well bound. Even though a removal of an imaginary part of the HF potential slightly affects self-consistency, 
solutions for unbound nuclei can easily be obtained. Other aspects of the GHF procedure remain unchanged 
as compared to the standard HF procedure for bound systems. 

Within the mean-field theory, pairing correlations can be included by means of the Hartree-Fock-Bogoliubov 
method (HFB). Its resonant-state extension is called Gamow-HFB |140j . For a bound nucleus, upper and lower 
components of one-quasi-particle Gamow states meet outgoing boundary conditions |140|, I141j . While the 
Gamow-HFB approach allows for the study of the structure of resonant states in weakly bound paired systems, 
its extension to unbound nuclei still needs to be worked out, even though an approximate scheme has been 
devised in Ref. [Hg. 

Another complex-energy, self-consistent approach that can be used to study particle unstable nuclei is the 
complex-scaled HF (CSHF) method |129j . In this approach, the HF potential is complex-rotated (j27p in order 
to have all its occupied states square-integrable. Contrary to GHF, self-consistency is fully maintained for 
unbound nuclei, because the CSHF potential is complex. However, one cannot extend CSHF beyond the HF 
level. For that, scattering states would have to be calculated in all directions of the complex energy plane 
in order to extend CS to non-resonant continuum. This is impossible because the CSHF potential diverges 
exponentially in some regions of the complex plane. Consequently, as CSHF cannot be used to generate a GSM 
basis, and pairing correlations cannot be treated within CS-HFB, the range of applications of CSHF is rather 
limited. 



4.5 Numerical tests of the Berggren completeness relation 

In this Section, we illustrate the validity and accuracy of Eq. ([6j) by means of several numerical examples. To 
this end, we consider the expansion of resonant states generated by a WS potential VE{r) (j3ip with a Berggren 
basis generated by a WS potential VB{r) of different Vq strength, but identical for all remaining parameters. The 

s.p. Hamiltonian to be diagonalized is represented by a continuous matrix in the eigenstates oihs = \-Vb- 

2m 

{un\- VVe\up) = -T-^Snp+ Un{r)A{r)up{r) dr, 

2m 2m Jo 

{u{k)\^ + VE\u{k')) = -—5{k-k')+ u{k,r)A{r)u{k' ,r) dr, 
2m 2m Jo 

{un\^ + VE\uik)) = / Unir)A{r)uik,r) dr, (48) 
2m Jo 

where Ui{r) and u{k,r) are, respectively, resonant and scattering states of the Berggren basis generated by /i^, 
and A(r) = VE{r) - VB{r). 

Following the contour discretization described in Sec. 14. 3^ the continuous representation of (j48j) is approxi- 
mated by a complex symmetric matrix: 

fp. h'^k'^ /•+00 

+ VE\ui>) = —^5ii> + / Ui{r)A{r)ui>{r) dr, (49) 

ZTTi Z771 Jn 
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where Ui{r) are discrete vectors representing both resonant and scattering states. The eigenvectors of ()49p 
represent the resonant and discretized scattering states of Ve- 

As an illustrative example, let us consider the expansion of a neutron 2^*3/2 eigenstate of a WS potential (j3ip 
in a Berggren basis. The radial wave function of the expanded state, denoted as UE{r), is a linear combination 
ofp3/2 basis states UB{k,r): 

u,(r)= E c,.uMK,r)+[ e(k)uMk,r),k, (50) 

ne{b,d) 

where Cn and c{k) are the expansion coefficients. As discussed in Sec. 14. 3| has to be discretized. 
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Figure 5: Distribution of the squared amplitudes c^{k) of the neutron state 2^3/2 of the WS in the basis 
generated by another WS potential. The amplitudes of both real (solid line) and imaginary (dashed line) parts 
of the wave function are plotted as a function of Re [A:]. The height of the arrow gives the squared amplitude 
of the 2p3/2 resonant state contained in the basis. The complex contour corresponds to three straight 
segments in the complex A:-plane, joining the points: A;o=0.0-iO.O, A;i =0.2-10.2, /c2=0.5-i0.0, and A;3=2.0-i0.0, all 
in units of fm~ . Top: Vo=59 MeV, =62 MeV. Here the expanded state is a 2^3/2 resonance. Bottom: 

Vo=64 MeV, Vq^^=59 MeV. Here the expanded state is bound. (Example taken from Ref. |119j where details 
of the calculation can be found.) 

The expansion coefficients (I50p are shown in Fig. [5] where the depth of the WS potential is varied to generate 
a resonant or bound 2^3/2 state. In both cases, the importance of the scattering continuum is noticeable. This 
example nicely demonstrates the invalidity of the pole approximation to provide a quantitative description of 
weakly bound or resonant states. 

More examples of s.p. expansions can be found in Refs. [119] (neutrons) and [53] (protons). Let us only 
mention that in the proton case the role of the scattering continuum is somehow reduced due to the confining 
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effect of the Coulomb barrier, especially at small energies. Still, for a precise description of eigenstates, the 
scattering component needs to be included. 

5 Gamow Shell Model 

5.1 Many-body GSM basis 

The discretized basis (ji2]l can be used to extend the completeness relation to the many-body case, in a full 
analogy with the standard SM in a complete discrete basis, e.g., the HO basis. The many-body basis of GSM 
corresponds to Slater determinants (SD) built from one-body states constituting the Berggren ensemble: 

\SD)n = \uij^---Ui^), (51) 

where \ui) are properly normalized one-body states (either resonant or scattering) of a generating potential. 
The completeness of the Berggren ensemble guarantees the closure relation for the many-body GSM basis: 

~1, (52) 

n 

where n runs over all possible SDs (I5ip . The approximate equality in ()52p is a consequence of the continuum 
discretization. 

5.2 Determination of many-body resonant states 

The GSM Hamiltonian matrix in the basis (j52p is complex symmetric. In a standard SM, the Lanczos method 
is often used to determine the eigenstates of the Hamiltonian matrix. It allows calculating selected states, e.g., 
ground state and low-lying excited states, in very large configuration spaces. However, the Lanczos method 
cannot be applied mutatis mutandis to the GSM matrix problem, due to the presence of scattering states. As 
many-body resonant states are embedded in the discretized continuum of many-body scattering states, the 
lowest-energy principle cannot be used to guide our choice. 

One way of selecting many-body resonant states was proposed in Refs. [1171 I143j dealing with systems 
having two valence nucleons. Here, the selection of two-particle resonances is done with a suitable choice of 
the contour L"*". It was realized that by choosing a square-shaped contour, the resonant states appear in the 
so-called allowed region of the complex-energy plane [11 7j and this facilitates their identification. 

In Fig. [6] the allowed region is a square-shaped area defined by a s.p. contour. Due to the low density of 
two-particle scattering states in this region, the trajectories of two-particle resonances can be easily followed 
as the strength of the residual interaction increases. It has been observed that the two-body resonant states 
are mostly based on the pole-pole (p-p) configuration, although the role of discretized continuum states (p-c, 
c-c) in Fig. [6] is not negligible. The result that the correlated two-particle resonant states are built mainly 
from configurations formed from s.p. resonant states gives strong support for the use of the overlap method 
discussed below. Unfortunately this convenient method might not be generalized easily to the case of more 
than two valence particles. 

Another solution to the problem of resonant-state selection, particularly useful in the case of several valence 
nucleons, has been proposed in Ref. |119j . It is based on the overlap method, utilizing the property that the main 
components of many-body resonant states are built from resonant SD. In the first step, the GSM Hamiltonian 
is diagonalized in a smaller basis consisting of s.p. resonant states only (a pole approximation). Here, some 
variant of the Lanczos method can be applied. The diagonalization yields the first-order approximation to 
many-body resonances |\I'o)- In the second step, |^'o) is used as a Lanczos pivot. Diagonalization of the GSM 
Hamiltonian in this subspace yields many-body eigenstates. The resonant state is the eigenstate |^) which 
maximizes the overlap |(^o|^)|- This method identifies the many-body resonant state provided its coupling to 
the non-resonant continuum is weak, i.e. 10-30%. This condition is usually satisfied, especially if the one-body 
basis has been optimized by way of the GHF procedure. 

Fig. El illustrates the overlap technique for J'^ = O"*" states in ^'^O considered as a four-neutron system outside 
the ^^O core [119j . Diagonalization of the GSM Hamiltonian matrix provides complex-energy eigenstates. 
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Figure 6: A typical square-shaped contour in the complex energy-plane used for calculation of two particle 
resonant states in Ref. |117j . Stars mark the s.p. pole (p) states, while the thick line denotes the scattering 
contour (c). Full dots denote the two-particle energies of the pole-pole type and those in which at least one 
single particle is in a scattering state. (See Ref. [ITTj for details.) 

Compared with Fig. [6l the number of many-body states is much larger and the regular pattern of non-resonant 
states reflecting the structure of the scattering contour is gone (the figure represents the projection of four- 
dimensional trajectories onto two dimensional space). While the two lowest (bound) states can be simply 
identified by inspection, for the higher-lying states it is practically impossible to separate the resonances from the 
non-resonant continuum. However, the procedure outlined above makes it possible to identify unambiguously 
the many-body resonance states. It is seen that the resonant states shown in Fig. [7] they all have negligible 
or very small positive widths, as expected. The many-body resonant states appear to be stable with respect 
to small changes of the contour (as the physical solutions should not be dependent on the deformation of the 
basis). On the other hand, the states representing the non-resonant many-body continuum move in the complex 
energy plane with the contour's deformation. 

5.3 GSM Hamiltonian 

Hamiltonian matrix elements between resonant states can be calculated with complex scaling (see Sec. [3]). 
However, this is not the case with scattering states, for which complex-scaled integrals may not converge 
regardless of the angle of rotation, as is the case for the normalization of scattering states (see Sec. 14. 2p . 
In the latter situation, since the result is analytical (Dirac delta), the problem can be dealt with by means 
of the contour discretization (see Sec. 14. 3p . However, it is very difficult to calculate matrix elements of 
general interactions; hence, in the first applications, only localized Hamiltonians were used. Examples are the 
Gaussian separable interaction [1171 1139j , the Surface Delta interaction [1181 1119j , and the Surface Gaussian 
interaction (SGI) [53j (see Sec. 16. ip . For these schematic forces, the calculation of two-body matrix elements 
is straightforward. In order to deal with realistic interactions, which are very complicated to express and 
infinite-ranged in shell-model coordinates, a new approach, based on a HO basis expansion technique, had to 
be devised (see Sec. 16. 2p . 

5.4 Matrix elements of electromagnetic operators 

In the GSM, the electromagnetic (EM) transition selection rules and the angular momentum and isospin 
algebra do not change compared to the standard SM. However, to calculate the EM transitions, one can no 
longer use the long wavelength approximation because of the presence of the non resonant continuum. Indeed, 
for the diagonal EM matrix elements {niiijiki\0\nfifjfkf) (ki = kf) between the scattering states, the complex 
scaling cannot be carried out. Furthermore, one cannot employ the long wavelength approximation to the 
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Figure 7: Complex energies of the O"*" states in ^''O calculated in the GSM+SDI model. One- (In) and two- 
neutron (2n) emission thresholds are indicated. The many-body resonant states (bound states and narrow 
resonances) are marked by squares. Note that all those state have zero or positive widths, as expected. The 
remaining eigenstates represent the non-resonant continuum. 

EM operators as they behave like r^. Without the long wavelength approximation, however, the continuum- 
continuum matrix elements become finite, because it is always possible to carry out a complex scaling with the 
Bessel function of the photon jL((7r), as g 7^ 0. As all the other matrix elements can be regularized, the EM 
matrix elements are all well defined. 

5.5 Optimization of the scattering contour of GSM using the Density Matrix Renormal- 
ization Group method 

The application of the Berggren ensemble to GSM is associated with the explosive growth in the number of 
configurations with the number of active particles and the size of the s.p. space. To ensure completeness of the 
Berggren basis, for each resonant s.p. state one should include a large number of s.p. scattering states 

{(i, j)^-^^; f = 1, • • • , M}c lying on a discretized contour . All these states become new active shells; hence, the 
configuration space of GSM grows extremely fast with the number of contour discretization points. Moreover, 
the resulting complex-symmetric GSM Hamiltonian matrix is significantly denser than the Hamiltonian matrix 
of a conventional SM. Hence, taming the dimensionality growth in solving the many-body Schrodinger equation 
is crucial. 

There are several specific features of the dimensional catastrophe in GSM that may facilitate this task. 
Firstly, of all GSM eigenstates, of major interest are the many-body resonant states. As discussed in Sec. 15.21 
these eigenstates represent a narrow class of spatially localized states, immersed in a sea of non-resonant 
scattering states. Eigenstates of the GSM Hamiltonian in the pole space can serve as the reference vectors. 
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helping to identify resonant eigenstates out of a huge space of ah GSM eigenvectors. Hence, instead of a direct 
diagonahzation of the GSM Hamiltonian matrix, a strategy should be developed for an optimization of the 
scattering space. 

The Density Matrix Renormalization Group (DMRG) method was introduced in Ref. |144j to overcome the 
limitations of Wilson-type renormalization groups to describe strongly correlated ID lattice systems with short- 
range interactions (see also recent reviews |1451 1146] ). Up to now, most of the DMRG studies were focused on 
equilibrium properties in strongly correlated closed quantum systems (CQS) with a Hermitian density matrix. 
Nevertheless, non-equilibrium systems involving non-Hermitian and non-symmetric density matrices can also 
be treated in the DMRG [147j . The DMRG approach is ideally suited to optimize the size of the scattering 
space in the GSM problem as the properties of the non-resonant shells vary smoothly along the scattering 
contour. The main idea of the GSM-I-DMRG truncation algorithm [1481 1149j is to gradually consider different 
s.p. shells of the discretized non-resonant continuum in the configuration space and retain only A'opt optimal 
states dictated by the eigenvalues of the density matrix with the largest modulus. Since the GSM Hamiltonian 
is complex symmetric, it is the complex (or generalized) variational principle [1501 115H I152j than the usual 
variational principle that governs the GSM-I-DMRG algorithm |149j . 

The iterative GSM-I-DMRG procedure of constructing the eigenvectors is divided into two phases: the warm- 
up phase and the sweeping phase. There are two essential conditions which guarantee that the GSM-I-DMRG 
method yields correct results. The first one is to ensure that all possible couplings in the many-body wave 
function, allowed by the symmetries of the problem, are present in the warm-up phase. This means that one 
should keep all the essential partial waves which are present in DMRG wave functions, even though they may 
be absent at the pole approximation level. The second condition is to ensure that the reference state 
which is a key ingredient in selecting the target state \^j) from the set of all GSM solutions of a given angular 
momentum and parity at each iteration step, is correctly described. 

The DMRG configuration space can be initially divided into two subspaces: A (the reference subspace built 
from all s.p. resonant shells {(?^l,^l, Ji), (n2,£2; J2); • • •} and one representative non-resonant shell for each 
remaining partial wave {(^p, Jp), {^r,jr), . . .} in the valence space) that is essential to produce many-body cou- 
plings, and B (the complementary subspace built from the s.p. non-resonant shells {{ii,ji)}c, {(^2, J2)}c, • • •)■ 
One begins in the warm-up phase by constructing all product states \k)A forming the reference subspace A. 

The many-body configurations in A can be classified in different families {n;jA} according to their number 
of nucleons n and total angular momentum j^. All possible matrix elements of sub-operators of the two-body 
Hamiltonian acting in A are calculated and stored. The Hamiltonian is then diagonalized in A to provide 
the zeroth-order approximation to the targeted eigenstate \J'^), called the reference state In the next 

step, the subspace of the first scattering shell belonging to the discretized contour L"*" is added. Within 

this shell, one constructs all possible many-body states, denoted as \i)B and grouped in {n^; j^} families, and 
all matrix elements of sub-operators of the Hamiltonian. This ensemble serves as a basis in which the GSM 
Hamiltonian is diagonalized. The target state \^j) is selected among the eigenstates of the Hamiltonian as the 
one having the largest overlap with the reference vector Based on the expansion: 

1*^) = E ^^a^WkAUA)) » \iB(jB))y , (53) 

by summing over the reference subspace A for a fixed value of Jb , one obtains the reduced density matrix |153j : 

P^^i'^{nB,3B) = 2.c.^o-5)C,^o-5) (54) 

kA 

which is complex-symmetric in the metric defining the Berggren ensemble. In the warm-up phase, the reference 
subspace becomes the 'medium' for the 'system' part in the B subspace. 

Truncation in the system part is dictated by the density matrix. The reduced density matrix is diagonalized 
and one retains at most A'opt eigenstates of p^: 

P^inBjB)\a)B = Wa\a)B (55) 
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having the largest nonzero values of \wa\- Then, one expresses the eigenstates \a)B in terms of the vectors \i)B 
in B: 

\a)B = Y.d?\i)B (56) 

i 

and all matrix elements of the sub-operators in these optimized states, 

B{a\0\P)B = E d?4 B{i\0\i')B (57) 

are recalculated and stored. In a similar way, the warm-up procedure continues by adding to the system part 
the configurations containing particles in the next scattering shells until the last shell in B is reached. At this 
point, all s.p. states have been considered, the warm-up phase ends, and the sweeping phase begins (cf. Refs. 
[14811149] for details). The sequence of sweep-up and sweep-down phases continues until convergence for target 
eigenvalues is achieved. 

It has been demonstrated [148^1149] that the fully converged GSM results, with respect to both the number 
of sweeps and the number of shells in the discretized continua, can be obtained using the GSM-I-DMRG 
algorithm. Several attractive features have been found. Firstly, the rank d'^'^^ of the biggest reduced density 
matrix to be diagonalized in GSM-I-DMRG is practically constant with respect to the number of shells A'sh- 
Hence, the ratio d^l^^/D, where D is the dimension of the GSM Hamiltonian matrix, decreases rapidly |148j . 
Secondly, the converged complex eigenenergies depend weakly on the number of vectors at each step of the 
iteration procedure, and the uncertainty of calculated eigenenergies weakly depends on N^h- This means that 
the convergence features of the GSM-I-DMRG procedure can be tested by varying Nsh and A'opt- Once those 
parameters are optimized in 'small-scale' GSM-I-DMRG calculations, the final calculations can be performed in 
the large model space to obtain fully converged results. Thirdly, the GSM-f-DMRG energy averaged over one 
sweep exhibits excellent exponential convergence with A'^opt, which allows to deduce the asymptotic value with 
good precision. These encouraging features of the GSM-I-DMRG algorithm open the possibility for systematic 
and high-precision studies of complex, weakly bound, or unbound nuclei, which require large configuration 
spaces for their correct description. 



6 GSM applications 

6.1 GSM applications with schematic Hamiltonians 

In the first applications of GSM, schematic core-|- valence Hamiltonians were applied. In Ref. |119j . neutron-rich 
isotopes of He and O have been studied by approximating the residual interaction in ()44l) with the Surface 
Delta Interaction (SDI): 

V{ri, Tj) = -VsDiS{ri - Ro)5{rj - i?o), (58) 

where rj = |rj|, Rq is the radius of the WS core potential, and Vsdi is the SDI strength. The main advantage 
of the zero-range SDI interaction is that its two-body matrix elements can be easily calculated and no complex 
scaling is involved. In the He case, the Berggren basis consisted of 0^3/2 and Opi/2 resonant states and the 
associated scattering continua and Pi/2- For the O isotopes, valence space consisted of the bound 0^5/2, 
lsi/2 states, the 0^3/2 resonance, and the ^3/2 scattering contour. Since other scattering channels can only 
impact real energies, they were not included as they could effectively be taken into account by renormalizing 
the SDI interaction strength. Results for the chain of oxygen isotopes are shown in Fig. [HI 

An early example of spectroscopic studies with GSM is shown in Fig. [9] which displays the level scheme of 
^^O calculated with GSM-I-SDI. The main experimental features are reproduced. In addition to energy levels, 
calculations were carried out for electromagnetic E2 and All transition probabilities. 

The GHF technique of basis optimization cannot be used effectively with SDI. Hence, in order to improve 
convergence of results, it was necessary to employ another interaction. For that, a delta component of the SDI 
interaction was replaced with a Gaussian form factor, thereby providing finite-range, while the second delta 
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Figure 8: Experimental (EXP) and predicted (GSM+V^^;/) |119j binding energies of 18-220. The energies are 
given with respect to the ^^O core. 



function remains in order to keep the interaction surface-peaked. The resulting surface Gaussian interaction 
(SGI) reads: 



d{r^ + rj-2Ro), (59) 



where is the range, V^qj is the strength in the JT channel, and all other values are defined similarly to 
Eq. (|58p. Two-body matrix elements can still be conveniently calculated with the SGI interaction, as its radial 
components are one-dimensional finite integrals. The fact that only one delta function appears in Eq. ()59p 
implies that the GHF-SGI potential is well defined in coordinate space, albeit of non-local character. Non- 
locality of spherical potentials is conveniently treated via the equivalent potential method [53], in which a 
non-local potential is replaced by a local but state-dependent potential. The results of GSM-SGI calculations 
for the neutron-rich He isotopes are shown in Fig. [TOj One obtains a reasonable description of the ground- 
state(g.s.) energies and excited states relative to the g.s. energy of ^He. 

A similar calculation for the lithium chain [53j suggests that monopole corrections have to be added to the 



(J T=0) 

T = matrix elements of SGI, modeled through a particle number dependence of VgQj- . While experimental 
data for ^Li and ''Li are fairly reproduced, the drip-line nuclei ^''Li and ^^Li are predicted to be overbound. 
This is due to a rather limited configuration space used in Ref. [53|; most importantly, the neglect of s-waves. 
Clearly, for a qualitative description of the Borromean halo nucleus ^^Li, an extended valence space and more 
realistic interactions are needed, as well as a proper treatment of the center-of-mass motion. 



6.2 GSM with realistic interactions 

The inclusion of realistic interactions in GSM has been a theoretical and practical challenge. Standard SM 
methods rely on the G-matrix formalism [154J and Lee-Suzuki similarity transformation |155j . The Lee-Suzuki 
projection operator formalism or the renormalization group technique can be used to derive low-momentum 
interactions, usually denoted as Viow-k [156J . The extension of the Lee-Suzuki formalism to the complex-energy 
GSM framework has been done in Ref. |139j where applications to Gaussian-type interactions can be found. 

Contrary to schematic interactions used in early GSM calculations, realistic interactions have a very com- 
plicated structure [157^ I158| I159j . They are defined in a center-of-mass frame, so that the transformation to 
the laboratory frame has to be carried out. For a HO basis, this can be done by using the Brody-Moshinsky 
coefficients |160j . However, the generalization of this scheme to arbitrary bases through the vector-bracket 
transformation [1611 1162^ 1163] is very time consuming, as the finite sums in Moshinsky transformation must 
be replaced by sums of two-dimensional momentum integrals. For the multipole decomposition of the realistic 
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Figure 9: The GSM level scheme of O calculated in the full sd space of Gamow states and employing the 
discretized (10 points) non-resonant continuum. The dashed lines indicate experimental and calculated 
one-neutron emission thresholds. As the number of states becomes large above the one-neutron emission 
threshold, only selected resonances are shown. The electromagnetic transitions discussed in |119| are indicated 
by arrows. (See Ref. |119j for details.) 

interaction, radial two-dimensional integrals have to be calculated with complex scaling. Unfortunately, for 
some integrals, CS is impossible to apply. 

The solution to all these problems has been proposed in Ref. [138j . It is based on the idea that the nuclear 
interaction V is localized in space around the nucleus, so that the HO expansion should quickly converge for 
bound states and narrow resonances. Expanded in a HO basis, two-body matrix elements of V can be written 
as: 

ax 

{ab\V\cd) = {aP\Vh6){a\a)m{ch){d\6), (60) 

where A^max is the number of HO shells, and Greek and Latin letters refer, respectively, to HO and Gamow 
one-body basis states. As V appears only in HO matrix elements (I60p . Moshinsky transformation can be 
easily performed. V can also be expressed in momentum representation in which HO states are analytical. 
Furthermore, Gamow states appear only in overlaps {a\a). The latter ones can be calculated without the need 
for complex scaling because of the Gaussian fall-off of HO states. 

The only approximation concerns the number Amax of HO shells used. Convergence for Amax — ^ +oo is in 
fact weak, i.e., while the energies and vectors of resonant states properly converge, the matrix elements (|60p do 
not. Note that this scheme is very different from HO expansion of SM states. In the standard SM, both kinetic 
and potential terms are expressed in a HO basis, whereas the method of Ref. |138j invokes HO expansion for 
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Figure 10: Experimental (EXP) and predicted (GSM) binding energies and spectra of helium isotopes obtained 
with the SGI Hamiltonian. The resonance widths are indicated by shading. The energies are given with respect 
to the core of ^He. The ground-state width of ^He is not known experimentally; hence, it is not indicated. 

nuclear interaction only. Since the kinetic and Coulomb terms are treated exactly in the asymptotic zone, the 
new method is capable of describing weakly bound and unbound states. 
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Table 2: Energy convergence for Oi^ and 2i^ states in He and 42 state in O versus the number A'^max of 
HO nodes in the expansion of the realistic low-momentum N'^LO nucleon-nucleon interaction. The imaginary 
part of the energy is shown in parentheses. A cut-off of 1.9 fm~^ was used for the center-of-mass momentum. 
The oscillator length is 6=2 fm. All energies are in MeV. (From Ref. |138j .) 

First applications of GSM with realistic interactions were carried out in Ref. |138j with the low-momentum 
potential Viow-k |156j for ^He and ^^O, viewed as two-neutron systems. The ^He or ^^O cores were approximated 
by a GHF potential. The same GHF potential was also used to generate the Berggren basis in spd neutron 
partial waves. Viow-k was obtained from the chiral N^LO interaction |159j with a momentum cut-off A = 1.9 
fm~^. The calculated energies of selected states in ^He and ^^O are shown in Table [6^ as a function of the 
HO states used in the expansion ([60|) . One can see that the results are already stabilized at A''max~10, even for 
resonant states having a sizeable width. In order to show that asymptotic properties can be reliably calculated 
in this scheme, Fig. [11] shows the density distribution of the valence neutrons in the halo nucleus ^He. The 
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results are already reliable at 



max — 



4 while perfect convergence is attained at A'^max=10. 












5 

r (fm) 



Figure 11: Neutron halo density of the ground state of ®He calculated in GSM+Viow-k [138j . The full 
convergence is achieved at A'max ~10, see the inset. 

In the context of ab initio many-body techniques, the Berggren ensemble has recently been employed in the 
coupled-cluster approach |164j to study neutron-rich He isotopes. This work demonstrates the usefulness of a 
complex-energy framework for a microscopic description of weakly bound and unbound systems. 

6.3 Center-of-mass treatment in GSM 

In SM approaches, special attention is paid to the removal of the spurious center-of-mass (CM) motion [165^166] . 
and GSM is no exception. The standard SM treatment of CM is based on the Lawson method [160j which, 
unfortunately, cannot be used within the GSM framework due to the fact that the matrix elements of (R 
is the CM coordinate) cannot be regularized using the complex scaling technique. One possible solution is 
to replace the CM Hamiltonian by the square of the component of the CM angular momentum operator 
L = R X P. Using of instead of L^^, would be equivalent but more complicated because is a four-body 
operator. In this generalized Lawson method, the CM part of SM nuclear states is not projected onto a Os HO 
state, but on a L= CM state. 

Another possibility is to solve the GSM problem directly in relative coordinates. While the Jacobi coordi- 
nates have to be excluded for practical reasons, as effects due to antisymmetrization become quickly intractable 
when the number of coordinates increases, one can define simple CM coordinates in a core+valence framework 
within the cluster-orbital shell model (COSM) [167j . In COSM, all coordinates are taken with respect to the 
core CM, so that translational invariance is strictly preserved. With COSM coordinates being particle coor- 
dinates. Slater determinants can be easily defined within this scheme. The translationally-invariant COSM 
Hamiltonian can be written as: 



where Ac and is the number of core and valence particles, respectively. The principal difference between the 
previously introduced GSM Hamiltonian (j44p and that of COSM is the term proportional to Pi • pj which takes 
into account the recoil of the active nucleons. The recoil term can be treated in the momentum representation 




(61) 
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using the Bessel expansion method and the same complex contour for all partial waves. Another possibility 
is to expand Pi in a basis of bound states, as done for realistic interactions in Sec. 16.21 

Because the transformation from the laboratory to COSM frame is very difficult to handle [168] . there are 
fundamental problems with using realistic interactions in a COSM framework. Nevertheless, this scheme is 
suitable for the development of effective interactions for light nuclei. 



6.4 Overlap integrals and spectroscopic factors in GSM 

Cross sections of direct reactions can be expressed as products of two factors. The first one depends on kinetic 
aspects of the reaction while the other one, depending on the structure of the states involved, is called the 
spectroscopic factor [169^ 1170^ I171j . The experimental deduction of spectroscopic factors can sometimes be 
difficult due to probe- and model-dependence [1721 11731 I174| . In this context, GSM is a useful tool as it 
provides both configuration mixing and a proper description of the particle continuum of both resonant and 
non-resonant character. 

Spectroscopic factors are defined as the norm of the one-nucleon radial overlap integral uej{r) \nO\ \171\ : 



u,,{r) = {^'/\ [{^"/S,') |£,i)J"'^), (62) 

where l^^'*) and are wave functions of nuclei A and A — 1 and \i,j) is the angular-spin part of the 

channel function. The angular-spin degrees of freedom are integrated out in Eq. (I62p so that uij depends only 
on the relative radial coordinate of the transferred particle. 

Using a decomposition of the channel in the complete Berggren basis B, one obtains: 

aUBmt\') {rij\B), (63) 



Ja 




S' =^{K-\\ai.{Bmit\'f, (64) 
where a'j'j{B) is a creation operator associated with a s.p. Berggren state \B). Since Eqs. (I63|64p involve 

£ ■ 

summation over all discrete Gamow states and integration over all scattering states along the contour , the 
final result is independent of the s.p. basis assumed. This is in contrast to standard SM calculations where the 
model dependence of SFs enters through the specific choice of a s.p. state a^^^- 

Of particular interest is a non-analytic behavior of direct reaction cross section near the particle emission 
threshold \175\ \176\ \177\ \178\ I179j . It manifests itself by a discontinuity in the cross section and/or of its 
derivatives (the so-called Wigner cusp); it often arises from a redistribution of particle flux when a new reaction 
channel opens. This phenomenon has been theoretically studied in Refs. \180\ \181\ \182\ \183\ \184\ \185\ [3T] and 
also recently in GSM [TMlfTST] . 

The spectroscopic factor pHe(g.s.)(8)p3/2]'^^ in ^He(g.s.), calculated in the GSM, is shown in Fig. [12] (top) as 
a function of the real energy ep^^^. The spectroscopic factor strongly depends on the position of the pole, 
especially around the In-emission threshold in ^He (6^3^2=0). As discussed in Refs. |186l 1187] , this dependence 
(in particular, a cusp due to a coupling with the ^He+n+n channel) follows the threshold behavior of the 
reaction cross section |175j . Specifically, the anomalous component of the spectroscopic factor below the In 
threshold in ^He behaves as (— e£j)^~^/^. Above the threshold, the spectroscopic factor is complex; the real part 
behaves as (e^j)^"*"^/^ while the imaginary part, associated with the decaying nature of ^He, varies as (e£j)^~^/^. 

Fig. [12] (bottom) illustrates the dependence of spectroscopic factors on the orbital angular momentum i of 
the transferred nucleon. Here, we consider the partial wave in ^^0 and the spectroscopic factor for the 
excited 0^ state of ^^O in the channel [^^0(3/2]^) (g> 4/2]°^ (see Ref. [W] for details). The behavior of S"^ is 

similar to that of (^IIe(g.s.)|[^He(g.s.) <8i P3/2]'^^)) except the variations are much weaker and the spectroscopic 
factor is continuous and smooth; it is its derivative that exhibits a cusp around the 0^3/2 threshold. This is 
again consistent with the general expectation that, for i=2, Re{S'^) below the In threshold of ^''O behaves as 
(— e^j)'^/^ while above the threshold Re(S'^) (/m(5^)) should behave as (e^j)'^/^ ((e^j)^/^). 
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Figure 12: Top: The real part of the overlap integral S'^ = (^He(g.s.)| pHe(g.s.)(E>p3/2]'^^)^ as a function of the 
energy of the 0^3/2 resonant state. ^He(g.s.) is bound for all values of eop.^^^. The solid line (GSM) shows the 
full GSM result. The dotted line (HO-SM) corresponds to the SM calculation in the oscillator basis of OP3/2 
and Opi/2 valence shells, and the dashed line (GSM-p) shows the GSM result in the pole approximation. The 
imaginary part of 5^ is shown in the inset. has been normalized so as to be equal to one for vanishing 
residual interaction. Bottom: Similar to the top portion except for the overlap in the excited Oj^ state of ^^O: 
(^^O(0;[)| [^^0(3/2+) (g) ^3/2]°"^)^ The first derivative of S"^ in the neighborhood of the 60^3/3=0 threshold is 
shown in inset (a) while inset (b) displays the imaginary part of S'^. 



To assess the role of the continuum, both resonant and non-resonant, the GSM results are compared in 
Fig. [12] with the standard SM calculations in a HO basis (HO-SM) and the results of a pole approximation 
(GSM-p). In contrast to GSM, spectroscopic factors in HO-SM and GSM-p vary little in the energy range 
considered, and no threshold effect is seen therein. A difference between GSM and GSM-p results illustrates 
the impact of the non-resonant continuum. Clearly, the continuum coupling should seriously be considered 
in studies of spectroscopic factors as the associated threshold effects can be as large as those due to, e.g., 
short-range correlations |188j . 

In Ref. |187j one-neutron radial overlap integrals were studied in GSM for various energy conditions of 
parent and daughter nuclei. It was concluded that Uij{r) can be very well approximated by a resonant state of 
a one-body potential which reproduces the complex generalized one-nucleon separation energy 

Sin{N) = E{N - 1) - E{N) = Sin{N) - [T{N - 1) - T{N)] . (65) 

Indeed, the asymptotic behavior of U£j{r) directly depends on Sin jl891 1190] : 

U£j{r) ~ e"*^^ , r — > +00 , k = \J 2mSin/h. (66) 

This is a straightforward generalization of a treatment of one-nucleon overlaps often used in SM studies where 
a real optical potential is employed with a depth adjusted to reproduce the experimental value of 5i„. 
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6.5 Comparison between GSM and COSM for light nuclei 



There have been quite a few apphcations of the Berggren ensemble to few-body systems. In the context of 
this review, of particular interest are COSM calculations of Refs. [191^ 1192] . They considered a complex-scaled 
Hamiltonian to study light drip- line nuclei such as ^~^He and ^^Li, described in terms of cluster+valence 
neutrons systems. The ^Li core of ^^Li is a complex nucleus which was described using a multi-configurational 
wave function [191]- The use of "V" and "T" types of coordinates in the three-body hybrid- TV model [193] 
("V" being related to COSM coordinates and "T" to cluster coordinates) makes it possible to take into account 
a large number of partial waves within a relatively small configuration space. 



COSM 



COSM 



1 



GSM 



(0^3/2)' 

(5'l)p3^2 

sum 

iOpi/2? 
('S'l)pi/2 

sum 



1.211-10.666 
-0.252+10.692 
-0.042-i0.026 
0.917 

1.447+10.007 
-2.658-10.042 
1.249+10.034 
0.038 



1.139-10.742 
-0.119+10.773 
-0.060-0.031 
0.960 

0.353-10.077 
-0.534+10.065 
0.221+10.012 
0.040 



1.105-10.832 
-0.060+10.881 
-0.097-10.050 
0.948 

0.226-10.161 
-0.198+10.224 
0.025-10.063 
0.053 



Table 3: Components of the ground-state wave function of ^He in GSM [139] and COSM [194J. 51 and S2 
indicate configurations with one and two particles in the non-resonant continuum, respectively. The COSM 
results obtained with partial waves up to ^max=5 are displayed in the second column, while the COSM results 
in a truncated space of and pi/2 partial waves are shown in the third column. (From Ref. [194j.) 



The COSM predictions for the He and O chains |194j have been compared to those of GSM [119\ 1139] . 
Table 16.51 shows the ground-state wave function decomposition for ^He obtained in both approaches. While 
configurations involving the 0^3/2 resonant state are similar, those associated with the Opi/2 broad resonant 
state are very different. This is due to the large coupling to the non-resonant continuum present in the pi/2 
partial wave with a weak contribution of the Opi/2 pole. While GSM is limited to and pi/2 partial waves, 
COSM includes all partial waves up to £=5, so that additional couplings change the relative occupancy of 
resonant and non-resonant pi/2 shells. This is nicely illustrated in a variant of COSM calculations in which a 
renormalized interaction was employed in a truncated space of partial waves with 1=1. As seen in Table [631 
GSM and COSM occupations are fairly close in this case [194] . 

One should stress, however, that the sum of components of the ground-state wave function of ^He over 
resonant state and non-resonant continuum states for p-^/2 and pi/2 partial waves is practically the same in 
both GSM and COSM. These quantities give the basis independent estimate of configuration mixing in the 
ground-state of ^He as given by GSM and COSM Hamiltonians, respectively. Hence, one can conclude that 
taking more partial waves does not change quantitatively the prediction about a structure of the ground-state 
wave function of ^He |119j . 

In many aspects, GSM and COSM are complementary. The use of complex-scaled Hamiltonian allows 
to calculate resonant states in COSM in a basis of discrete bound states. Exact asymptotic behavior, which 
is efficiently treated by means of hyperspherical coordinates, has to be ensured. In practice, however, this 
requirement is impossible to fulfill with more than three clusters. This seems to be a major limitation of 
COSM. The use of overlap method in GSM circumvents this problem, but the price to pay is blurring different 
emission channels, so that one cannot easily identify different partial widths. This problem is currently under 
investigation. 
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7 Perspectives 



Shell structure is a fundamental property of finite Fermi systems. In nuclear physics, the nucleonic shell 
structure is characterized by the appearance of large energy gaps in the s.p. spectrum. In open shell nuclei, SM 
explains the shell evolution in terms of the configuration mixing due to a residual interaction. Experimental 
discoveries over the last decade, accompanied by theoretical studies, have demonstrated the fragility of the SM 
concept, in particular the separation of the Hamiltonian into a s.p. potential (that generates shell structure 
with immutable magic numbers) and a two-body residual interaction [1]. As a result of configuration mixing, 
new magic numbers in neutron-rich nuclei appear far from stability, while some traditional ones disappear. 
Another systematic effect of large neutron excess, illustrated in Fig. [H is a decrease of one-neutron separation 
energy when moving towards the neutron drip line. 

GSM is the OQS formulation of SM for a self-adjoint Hamilton operator. This choice offers a number of 
conceptual advantages. For instance, the transition from a bound to unbound regime, either within a single 
nucleus or in the long chain of isotopes (isotones), can be viewed as an opening in the configuration space and 
described without changing the Hamiltonian. The many-body nuclear Hamiltonian in this formulation does 
not describe just one isolated nucleus {N,Z), but all nuclei and all nuclear states that are coupled through 
various decays and captures. This idealization offers a right physical picture of the many-body system and its 
interactions. 

Already early applications of GSM to the helium chain revealed the essential role of many-body contin- 
uum states in the binding mechanism of ^He (bound), ^He (unbound), and ^He (bound). In addition, they 
demonstrated that the configuration mixing involving continuum space may significantly affect occupancies of 
different s.p. orbits in many-body states in the vicinity of particle-emission thresholds, thus inducing specific 
variations of gaps in the spectrum of effective s.p. energies. The resulting picture differs even further from 
that of the standard SM: the s.p. shell structure probed by GSM becomes a fragile concept in weakly bound 
nuclei [1]. Studies of many-body wave functions of GSM inform us how subspaces of discrete and scattering 
states are intertwined. Having a complete basis which allows to describe bound, weakly-bound and unbound 
states on the same footing is the only way to guarantee the unitarity. This fundamental requirement is not 
respected in any CQS formulation of the many-body theory. As shown in several examples, restoration of 
unitarity can strongly affect both values and behavior of spectroscopic factors. The unitarity lies also at the 
heart of the Wigner threshold effect in the multichannel case. Further consequences of the unitarity restoration 
on low-energy reactions, mirror symmetry-breaking, pairing correlations, and multi-nucleon decays, deserve 
further investigations. 

Over the last few years, the Berggren ensemble and GSM became an accepted tool for nuclear structure 
studies. An efficient algorithm for finding discrete resonance states in the non-resonant continuum has been 
proposed |119j and a new method, based on the application of the DMRG technique, was developed to tame 
the explosive growth of the Fock space [148] . Recently, the Berggren ensemble has been applied in an ab initio 
description of open quantum systems in a complex-energy coupled-cluster approach |164j . The GSM serves 
also as a benchmark model to test other approaches to weakly bound states, such as the COSM |194j . Progress 
has been made in applying Berggren ensembles to the mean- field description of nuclei through Gamow-HFB. 

Many challenging problems remain that could be illuminated with GSM — thanks to its ability to follow a 
quantum transition from a bound-state to unbound regime. For instance, GSM is an excellent tool to explain an 
'alignment' of a many-body wave function with a channel wave function in the vicinity of a channel threshold 
|195j . This mechanism, which is a likely source of clustering effects seen near different threshold openings, 
can be studied in GSM. A near-threshold behavior of pairing correlations is an example from another domain 
where systematic GSM studies could be most helpful. That is probably a most transparent illustration of a 
near-threshold clustering. 

To improve the detailed description of bound and resonance states spectra in GSM, a development of 
effective interactions in RHS is an urgent need. This will open a possibility for investigations of configuration- 
mixing effects involving continuum states in certain isotope chains, e.g., ^Li (bound), ^'^Li (unbound), ^^Li 
(bound), and ^^Li (unbound). In this example, GSM could provide a quantitative estimate of admixtures 
involving coupling of two neutrons (one neutron) with different states of ^Li (^'^Li) in the ground state of ^^Li. 
Recent experimental studies of the neutron transfer reaction p(^^Li,^Li)t |196j provide a strong motivation. 
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Much effort has been devoted to studies of spectral degeneracies associated with the avoided level crossings, 
focusing mainly on the topological structure of the Hilbert space and the geometric phases [1971 1198] . Among 
these degeneracies, one finds exceptional points [1991 12001 1201j where two eigenfunctions and associated eigen- 
values become equal. Up to now, exceptional points have not been studied within RHS. In this respect, GSM 
could provide a first realistic and consistent many-body framework in which properties of wave functions close 
to a RHS singularity could be described. 

In summary, a complex-energy shell model opens a new window for unification of structure and reaction 
aspects of weakly bound or unbound nuclear states, based on the open quantum system framework. The 
developments reviewed here are not solely limited to nuclei; they can be of interest in the context of other open 
quantum systems in which the coupling to the scattering continuum is present. To describe such systems, one 
has to give up either the concept of the Hilbert space or the self-adjoint nature of the Hamiltonian if one wants 
to keep the simplicity and conceptual and practical advantages of the standard SM framework. 
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